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Abstract 

A generalization of Faddeev's approach of the 3-body problem to the many- 
body problem leads to the method of increments. This method was recently 
applied to account for the ground state properties of Hubbard-Peierls chains 
(JETP Letters 67 (1998) 1052). Here we generalize this approach to two- 
dimensional square lattices and explicitely treat the incremental expansion up 
to third order. Comparing our numerical results with various other approaches 
(Monte Carlo, cumulant approaches) we show that incremental expansions 
are very efficient because good accuracy with those approaches is achieved 
treating lattice segments composed of 8 sites only. 
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The understanding of properties of strongly interacting fermions has been an intense 
topic of research for the past decade, in part due to the interest in the properties of high 
temperature superconductors. As exact solutions are known only for selected integrable 
models, numerical methods gained importance to provide benchmarks for analytical ap- 
proaches which necessarily use approximations of all kinds. 

Incremental expansions have been used in quantum chemistry to account for properties 
of molecules and solids fl|]. In a recent work || these methods were combined with cumu- 
lant expansions to provide a solid footing for numerical implementations. One result was, 
that incremental expansions can be interpreted as approximative ways to solve Faddeev-like 
equations for the N-body problem. 

To explain this in more detail, let us consider a Hamiltonian 

H = H + J2H n ■ (1) 

n 

Assume that we can find the eigenenergies and eigenvectors of Ho and of Ho + H n for any 
n. Suppose further that S n is the scattering operator associated with H + H n . Then 
Faddeev-like equations read as M 

T n = S n (l + T,) , S = Y,T n , (2) 

where S is the scattering operator of the full Hamiltonian H . Assuming the S n to be small, 
one obtains in lowest order T n = S n and S = J2n S n . Having the S operator one can calculate 
e.g. the ground state energy of H (for details see @). 

The ground state energy E in zero order of the incremental expansion equals to the 
ground state energy E^ of Ho- In first order we take into account the effect coming from 
adding one single H n . Denote the ground state energy of H + H n by E n . Then the first 
order increment to the ground state energy is given by 1^ = E n — E^ \ This increment 
measures the change of adding H n to the ground state energy. Summing up all increments 
of first order (for different n) we obtain the ground state energy E^ to first order 

£(1)=£(0) . ( 3 ) 
n 
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Let us now calculate the second order increments, which are taking into account the simulta- 
neous effect of adding two terms H n and H m . Here one has to subtract the corrections coming 
from adding both terms separately. Denote the ground state energy of (H + H n + H m ) with 
E nm . Then the second order increment is given by ijj^ = E nm — — — . The 
ground state energy to second order E^ is then given by (see Chapter 5.2.3 in ||T|) 

^) =E (0) +E/ (D +E/ ( 2r ) . (4) 

n n<m 

This procedure can be easily continued to higher orders, and the ground state energy will be 
a sum over increments of all orders. It requires the exact calculation of E nm when going to 
second order ( respectively E nm i when going to third order etc). Clearly this procedure does 
not imply certain topological structures induced by J2 n H n , so that we are not restricted 
to certain space dimensions if considering spatially extended systems. Note that there are 
many ways to split a full Hamiltonian H into a trivially solvable part Hq and the terms H n . 
Each of these ways would generate its own incremental expansions. Finally these expansions 
are not restricted to the calculation of the ground state energy only, but can be also applied 
to other ground state properties and to excited states 0. 

This method was successfully applied to one-dimensional lattices [J3J. There the special 
topology of a Id system lead to the successive cancellation of lower order increments when 
proceeding to higher orders. Here we use this method to account for the ground state energy 
of the two-dimensional Hubbard model on a square lattice at half filling. The Hamiltonian 
in dimensionless units is given by 

H el = ( C lj,aCi+l,j,a + 4^*4+1,* + h - C ) + U n i,jA U i,U ■ ( 5 ) 

Here % and j are integers (which denote the x and y coordinates of the lattice points). To 
proceed we have to define a splitting of H e i into Hq and J2 n H n - Here we consider Hq as 
given by H e \ when all vertical bonds and each second horizontal bond are missing (see also 
Fig.l): 

H = ~Y1 ( C 2i+j,j,aC2i+j+l,j,a + h.C.) +Uj^ ^71^,1 . (6) 



Thus Ho is a set of horizontally alligned noninteracting dimers. Note that Hq already 
contains all correlation terms of (ffy. The terms H n are then the missing bonds i.e. all vertical 
bonds Y,i,j( c i,j,(T C i,j+i,(T J r h.c.) and the missing horizontal bonds J2i,j( c 2i+j-i,j,a c 2i+j,j,a + h.c.) . 

For sake of concreteness let us assume that our initial model Hamiltonian (|5|) has an even 
and finite number of 2N sites, and periodic boundary conditions. In zero order the ground 
state energy E of the whole system and the energy per site e are given by 

E® = NE(C ) , e(°) = l -E{C Q ) , (7) 

where E(Cq) is the ground state energy of the Oth configuration Co which per definition is 
a dimer with two electrons (Fig.2). 

The first order increment is given by adding one of the missing bonds. The ground state 
energy of this case is given by sum over the ground state energies of iV — 2 dimers (Co) and 
of an open segment of two coupled dimers C\ (see Fig.2): 

1(d) = E(d) + (N — 2)E(C Q ) - NE(C ) = E{C,) - 2E(C ) . (8) 

We will encode all configurations of linked dimers to be considered by C n . Note that all 
increments are independent of the given position of the returned bond because we assumed 
periodic boundary conditions. Next we need to account for the weight factor w(C\) of I(Ci), 
i.e. the number of increments per dimer of H having the same energy. It is easy to see that 
w(C\) = 3 (all possible realizations are shown in Fig.2). The ground state energy in first 
order is then 

E (i) = E (o) + 3NI{Ci) = N ( 3E ( Cl ) - 5£(C )) , e« = \{3E{Ci) ~ 5E(C )) . (9) 

Already at this stage, though the considered configurations are equivalent to those of a Id 
chain at the same order of incremental expansions we find a difference in the energy per 
site due to the increased number of nearest neighbours of the 2d lattice as compared to the 
Id case. 
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In second order we have to add two of the missing bonds. Nonzero contributions come 
from cases when the two returned bonds are linked 0. Then we have two nonzero config- 
urations C 2 (open chain with six sites) and C 3 (which is already incorporating topological 
effects of the 2d system) in this order (see Fig. 3). Their weight factors are 10(62) = 9 and 
w(C 3 ) = 6. The corresponding increments are 

I(C 2 ) = E(C 2 ) - 3E(d) - 21(d) = E(C 2 ) - 2E(C 1 ) + E(C ) , (10) 
I(C 3 ) = E(C 3 ) - 3E(C ) - 21(d) = E(C 3 ) - 2E(C l ) + E(C ) . (11) 

The ground state energy can be evaluated according to @: 

= ^E(C ) + h(d) + h(C 2 ) + 6I(C 3 ) = 5E(C ) - 13.5£(d) + A.hE(C 2 ) + 3E(C 3 ) . 

(12) 

Since the topology of the configurations starts to be different from those appearing in a Id 
system no trivial cancellation of lower order increments takes place anymore. 

In third order we add three missing bonds (again only configurations when all 
three returned bonds are linked do contribute). We obtain six different configurations 
C4, C 5 , C 6 , C 7 , C 8 , C 9 with corresponding weight factors 2, 27, 2, 10, 18, 32. They are shown 
in Figs. 3,4. The corresponding increments are 



I(C 4 ) - 


= E(C 4 ) - 


- I(C 2 ) - 


2/(C 3 ) 


- 31(d) 


-3E(C ) , 


(13) 




m) = 


E(C 5 ) - 


2I(C7 2 ) 


- 31(d) 


-4£(C ) , 


(14) 




m) = 


E(C 6 ) - 


3/(C 3 ) 


- 31(d) 


-AE(C ) , 


(15) 




I(C r ) = 


E(C 7 ) - 


2/(C 3 ) 


- 31(d) 


-4£?(C ) , 


(16) 


ncs) 


= E(C 8 ) 


- /(C7s) - 


- 2/(c 2 ; 


) - 3I(d] 


) - 4E(Co) 


(17) 




= E(C 9 ) 


- I(C 2 ) 




- 31(d) 


-4E(C ) . 


(18) 



It is too tedious to explicitely write down the formula for the ground state energy. In the 
following we will present the results of numerical calculations. 
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We use a Lancosz algorithm to compute the ground state energies of our considered 
configurations. In Fig. 5 we show the dependence of the ground state energy per site e 
on the interaction parameter U in third order of the incremental expansion. In the inset 
of Fig. 5 we show the [/-dependence of the different incremental contributions. We find 
that the contributions coming from 2nd and 3rd order are small compared to the Oth 
and 1st order. For free electrons U = we compare e for Oth, 1st, 2nd and 3d orders - 
-1; -1.708204; -1.768700; -1.6335775 with the exact number e(U = 0) = -1.621139 (see 
also Table 1 and Fig.5). This gives a relative error of only 0.8% ! For U = 1, 2, 4 we compare 
e in 3d order in Table 1 with quantum Monte Carlo (QMC) calculations of Moreo et al [|J, 
where lattices with sizes up to 16 x 16 were used and extrapolations were carried out. The 
relative difference (Table 1) is less than 2%. Note that the slight increase in the relative 
error with increasing U at least partially has to be attributed to the circumstance that the 
QMC calculations become less exact with increasing U. Also in Fig.5 the results of Polatsek 
and Becker on related projection operator techniques using cumulants are shown ||. These 
calculations take spin flips up to second order and charge fluctuations up to second order 
(in terms of our notations) into account. The agreement is very good in a broad range of U 
values. The exception is the limit of small U, where the projection technique becomes less 
accurate. Finally for large U the Hubbard model transforms into the Heisenberg model. Do- 
ing the same calculations for the Heisenberg model with J = 1 we compare e for Oth, 1st, 2nd 
and 3d orders - -0.375; -0.5490375; -0.578994; -0.6695330 with the results of QMC calcu- 
lations of the ground state energy of the Heisenberg model e ~ —0.669 by Runge (see 
also 0, 0, || and with similar results). The relative difference is less than 0.08% (see 
also Table 1). To conclude this part, we emphasize that our results yield a high precision 
(relative differences of the order of 2% and less) in the whole U range, which has not been 
achieved by any of the other methods discussed. 

Let us emphasize that the presented method is not just a clever way of making finite 
size extrapolations. To show that, we consider the two-dimensional antiferromagnetic Ising 
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model H = J2ijS-Sj on a square lattice with spin 1/2. The ground state energy per 
site is given by -1/2, which is simply the result of each site having two bonds, each bond 
contributing with an energy of —1/4. Any finite size calculation of this energy would deviate 
from the exact value, because it would involve the energy of spins at the boundary of the 
finite size cluster, where the number of contributing bonds per site is less than 2. However 
it is an easy task to check, that our method gives E{Cq) = —1/4, I{C\) = —1/4, I(C n ) = 
for 0,1. Thus our expansion terminates after the first order, and in this order we obtain 
precisely e = —1/2. This should make clear that incremental expansions are at any stage 
yielding results for the infinite lattice. 

Using exact diagonalizations we could extend the calculations even further up to 6th or- 
der, i.e. up to adding six missing bonds. This needs a careful classification of all contributing 
configurations and their weight factors. We are currently working on this project. Notice 
the extremely high precision which we achieve already in 3d order, where the largest systems 
we have to deal with consist of 8 sites. Without any extrapolation we obtain a precision 
which e.g. in QMC is achieved by considering systems with size up to 16 x 16 = 256 sites 
and additional extrapolations (cf. Fig. 5 in [ID]). The reason for this is the fact that we use 
a scheme which at each level describes an infinite system, and accounts for the important 
topological structures through the weight factors. This appears to be much better than just 
to consider a finite lattice with a certain size. Having the ground state energy with that 
accuracy, we plan next to account for the dimerization of a two-dimensional Hubbard- Peierls 
system. Work is in progress. 



We thank P. Fulde for continuous support and K. W. Becker, S. Denisov, P. Fulde and 
R. Hetzel for helpful discussions. 



7 



REFERENCES 

[1] P. Fulde. Electron Correlations in Molecules and Solids. Third enlarged edition. Springer 
Berlin, 1995. 

[2] K. Kladko and P. Fulde. Int. J. Quantum Chem., 66:377, 1998. 

[3] J. Malek, K. Kladko, and S. Flach. JETP Letters, 67:1052, 1998 (Pis'ma Zh. Eksp. 
Teor. Fiz., 67:994, 1998). 

[4] A. Moreo, D. J. Scalapino, R. L. Sugar, S. R. White, and N. E. Bickers. Phys. Rev. B, 
41:2313, 1990. 

[5] G. Polatsek and K. W. Becker. Phys. Rev. B, 54:1637, 1996. 
[6] K. J. Runge. Phys. Rev. B, 45:7229, 1992. 

[7] Z. Weihong, J. Oitmaa, and C. J. Hammer. Phys. Rev. B, 52:10278, 1995. 

[8] E. Manousakis. Rev. Mod. Phys., 63:1, 1991. 

[9] K. W. Becker, H. Won, and P. Fulde. Z. Phys. B, 75:335, 1989. 

[10] S. R. White, D. J. Scalapino, R. L. Sugar, E. Y. Loh, J. E. Gubernatis, and R. T. 
Scalettar. Phys. Rev. B, 40:506, 1989. 



8 



TABLE 1. 

Comparison of e from 3rd order increments with QMC results and exact value for U = 0. 



u 


e (ref) 


e (3rd order inc.) 





-1.621 (exact) 


-1.634 


1 


-1.376 ( §) 


-1.400 


2 


-1.172 ( §) 


-1.191 


4 


-0.841 ( 1) 


-0.856 


oo (Heisenberg) 


-0.669 ( §) 


-0.670 
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FIGURE CAPTIONS 



Fig.l 

Schematic representation of H . Circles - lattice sites, lines - bonds kept from H. 
Fig.2 

Different configurations (see text). 
Fig.3 

Same as Fig.2 
Fig.4 

Same as Fig.2 
Fig.5 

e versus U. Solid line - result from third incremental order; filled circle - exact value at 
U = 0; open squares - QMC results (4j; filled diamonds - projection operator results ||. 
Inset: Dependence of incremental contributions to e on U. Solid line - 0th order; thick 
dashed line - 1st order; thick long dashed line - 2nd order; dotted line - 3rd order. 
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Fig.1 



